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We calculate the zero-temperature gap and quasiparticle weight of the half-filled Hubbard model 
with a random dispersion relation. After extrapolation to the thermodynamic limit, we obtain 
reliable bounds on these quantities for the Hubbard model in infinite dimensions. Our data indicate 
that the Mott-Hubbard transition is continuous, i.e., that the quasiparticle weight becomes zero at 
the same critical interaction strength at which the gap opens. 



PACS numbers: 71.10.Fd, 71.30.+h 

As pointed out by Mott Q many years ago, a suf- 
ficiently strong electron-electron interaction produces a 
gap for current-carrying charge excitations in systems 
with an integer average number of electrons per lattice 
site. In their simplest form, his ideas can be understood 
within the one-band Hubbard model Q , in which the ki- 
netic energy of the electrons competes with their purely 
local interaction. When there is one electron per lattice 
site on average and the interaction strength, U > 0, is 
large compared to the bandwidth, W, the charge gap is 
of the order of A(U ^ W) = U — W, irrespective of a 
possible symmetry breaking in the ground state. Hence, 
in the absence of nesting, the Hubbard model at half 
band-filling contains a zero-temperature quantum phase 
transition from a metal to the Mott-Hubbard insulator 
at some critical interaction strength U c of the order of 
the bandwidth, U c w W §. 

The Mott-Hubbard transition poses an intricate many- 
body problem, and few exact results are available. In 
one dimension, the quantum phase transition can be in- 
vestigated in the Hubbard model with a linear dispersion 
relation in which the absence of nesting precludes an an- 
tiferromagnetic instability at U = + . In this model, the 
metal ceases to exist when the gap opens at U c = W j| , 
i.e., the Mott transition is continuous 

The Hubbard model cannot be solved exactly in d > 1 
dimensions. However, many simplifications arise in the 
limit d — *■ oo ||. In this limit, the Hubbard model can 
be mapped onto a single-impurity problem in a bath 
of electrons whose properties must be determined self- 
consistently (dynamical mean-field theory); for recent re- 
views, see Refs. From this mapping it can be 
shown that a metallic ground state must be a Fermi liq- 
uid whose quasiparticle weight Z is identical to the size of 
the jump discontinuity of the momentum distribution ||. 

Within the dynamical mean-field theory, the interact- 
ing Green function must be calculated explicitly at all 
frequencies in order to iterate the dynamical mean-field 
equations. In this work we avoid this self-consistency 



scheme altogether by studying the Hubbard model with 
a random dispersion relation, which is equivalent to the 
Hubbard model in infinite dimensions We work di- 
rectly with static ground-state properties such as the gap 
and the quasiparticle weight (see Eqs. (0,|^) below) which 
can be used to determine the nature of the Mott-Hubbard 
metal-insulator transition. We calculate these quanti- 
ties using exact diagonalization on finite lattices of up 
to L = 14 sites, and extrapolate to the thermodynamic 
limit. As shown in Fig. |l|, our extrapolated results for 
the gap and the quasiparticle weight are consistent with 
a continuous Mott-Hubbard transition at U c as W for a 
semielliptic density of states. 
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FIG. 1. The quasiparticle weight Z(U) and the 

one-particle gap, A(U), extrapolated to the thermodynamic 
limit separately for even and odd system sizes. 

Our results are in strong disagreement with the sce- 
nario of a preformed gap proposed by Georges, Kotliar 
et al. [Q. They find a metallic solution of the self- 
consistency equations up to some J7 C ,2, and an insulating 
solution down to £/ c .i < [4,2, with U c< 2 — U Ct i « 0.4VF. 
They argue that the transition occurs at U c = U c ,2 be- 
cause the metallic solution has lower energy. The quasi- 
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particle weight then vanishes continuously, but the gap 
opens abruptly fjj . A coexistence regime at zero temper- 
ature implies a line of first-order metal-insulator phase 
transitions at finite temperatures, and it is tempting to 
explain the qualitative features of the phase diagram 
of V2O3 solely on the basis of electron-electron interac- 
tions 0. Recently, however, the arguments presented in 
favor of a discontinuous transition have been scrutinized 
and argued to be insufficient Pj4|p0|] . 

The Hubbard model on an L-site lattice is given by 

L L 

H= ^2 til - m)c+ a c m , a + U^hi A ni A , (1) 

l : m— l;cr I — 1 

where o[ a (c; lCT ) creates (annihilates) an electron with 
spin a on lattice site I, h^ a — cf^ci^ denotes the num- 
ber operator on site and t(r) — t*(—r). The ki- 
netic energy operator is diagonal in momentum space, 
T = J2k,a £ kci irT c k , < y, e k = (l/L)J2 r t(r)exp(-ikr). 
Here, c\ a — \/i/LJ2i exp(ifc^)c i f CT creates an electron 
with Bloch momentum k. 

Two important functions which characterize the dis- 
persion relation are the non-interacting density of states, 
D(e) = (l/L)^2 k S(e — e k ), and the joint density of 
states, D q (e ll e 2 ) = (l/£)Efc^( £ i ~ £k)5{e 2 - e k+q )- 
The latter quantity gives the momentum-averaged prob- 
ability for the event that a momentum transfer q gen- 
erates a transfer of kinetic energy from E\ to e 2 . In 
the random dispersion approximation (RDA) Q the dis- 
persion e k is replaced by the random dispersion £ RDA 
with the same density of states, D RDA (e) = D(e), 
but a completely uncorrelated joint density of states, 
fi A (£i,£2) = D(ei)D(e 2 )- All higher correlation func- 
tions must also factorize correspondingly. In order to 
satisfy these conditions, we randomly assign an energy 
e fc DA chosen from the probability distribution D(e) to 
each momentum k. In this version of the RDA all forms 
of magnetic order are suppressed ("fully frustrated hop- 
ping" 0) such that magnetic order does not conceal the 
Mott-Hubbard transition. 

One way to see the equivalence of the RDA to the 
limit of infinite lattice dimensions is the diagrammatic 
approach. The value of the local non-interacting Green 
function in the RDA is of order unity, and identical to 
that of the Hubbard model in infinite dimensions for a 
given D(e). For all other lattice distances the Green 
function behaves like \fTJ~L Two vertices I and 

?7i of a skeleton diagram in the proper self-energy are 
connected by three independent Green function lines |^] 
such that the contribution from I — m ^ vanishes like 
LO{L- 3 / 2 ) = 0{L~ x l 2 ). Thus, the self-energy becomes 
purely local. Moreover, the skeleton expansion of the 
Hubbard model in infinite dimensions and in the RDA 
are identical. The local interacting Green functions also 
agree to all orders in perturbation theory since they de- 



pend on momentum only through e k . Thus, the self- 
energies of both approaches are the same, i.e., they lead 
to the same mean-field theory. 

For a numerical analysis we must specify the disper- 
sion relation of a realization Q of the randomness, 
which then defines the corresponding Hubbard Hamilto- 

nian H Q = J2 k ,a £ k^k,a^,c + U Y^i ^,T™U- Eacn H Q is 
a well-defined Hamiltonian on a formally one-dimensional 
lattice which can be treated using numerical techniques. 
Here we use exact diagonalization with the Davidson al- 
gorithm p"l[ to obtain energies and equal-time matrix 
elements of low-lying states which are accurate to close 
to machine precision. Since we will have to average 
over configurations Q which correspond to diagonal disor- 
der in momentum space, we perform the diagonalization 
in the momentum-space basis. We adopt antipcriodic 
boundary conditions for an even number of sites L and 
periodic boundary conditions for odd L. 

For a given U, L, we calculate the single-particle gap 

AS = E*{L + 1) + E$(L 1) - 2E${L) , (2) 

and the momentum distribution 

^ S = ^E(*o S l<^|*o Q >, (3) 

a 

for each configuration Q, where l^^) is the normalized, 
paramagnetic ground-state wave function at half band- 
filling, and E®(N) is the Af-particle ground-state energy. 
Since the self-energy is fc-independent, the momentum 
distribution depends on k only through Ek' nj? = n^(e k ). 
Thus, Z Q = n Q (s = 0~) - n Q (e = 0+) defines the quasi- 
particle weight for configuration Q at half band-filling. 

To obtain the mean value of a physical quantity, we av- 
erage over Rl selected configurations Q, e.g., A(U; L) = 
(1/Rl) ^2q A s ([7; L). Random sampling determines the 
mean value with accuracy SA(U;L) — <ta{U; L)/^/R l , 
where [a A (U;L)f = (1/ Rl)J2q[^ q {U ; L) — A(U\ L)] 2 
denotes the width of the distribution for given U and L. 
Finally, we must extrapolate the result to infinite system 
sizes, e.g., A(Z7) = lim^oo A(U ; L). 

Here we present results for a semielliptic density of 
states, D(e) = &/(irW)[l - {2e/W) 2 ] 1 ' 2 for |e| < W/2. 
In one dimension this density of states translates into 
a dispersion relation which obeys the implicit relation 
jfe/2 = (2e k /W) [1 - (2e k /W) 2 } 1 / 2 + arcsin (2e k /W). In 
the following we measure energies in units of t = W/(2ir). 

We choose a particular realization for a configura- 
tion Q as a permutation of the L energy values {s k } over 
the momenta {k}. This choice avoids fluctuations in the 
bare density of states. We then select those permuta- 
tions Q for which the joint density of states D^(ei,s 2 ) 
(approximately) factorizes. For small finite systems, 
we analyze the distribution of hopping matrix elements 
\tQ(r)\ 2 = (1/L) E 9 exp(z 9 r) £ fc e§e§ ±q . In the thermo- 
dynamic limit we find |t RDA (r)| 2 = t 2 = (l/£)£ fc 4, 
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independent of r. We exclude configurations Q which 
strongly deviate from this relation, retaining about ev- 
ery second (randomly chosen) configuration. For L = 
4,5,6, we calculate all unique configurations; we select 
R L = 500 for L = 7,8,9, R L = 200 for L = 10,11,12, 
and Rl — 50 for L — 13, 14. In order to improve the 
statistics for small system sizes, we choose independent 
random dispersions for spin-"f and spin-| electrons. We 
emphasize that the L — > oo extrapolation does not de- 
pend strongly on how the configurations Rl are chosen; 
our selection criteria simply reduce the statistical errors 
in order to make the extrapolation more accurate. 
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FIG. 2. One-particle gap A(i7; L) plotted as a function of 
1/L for various U values. Separate extrapolations are done 
for even L (solid lines) and odd L (dotted lines). 

In Fig. H we plot the one-particle gap A(U ; L), Eq. (||), 
as a function of 1/L for various values of U. The lines 
are separate 0(1/ L 2 ) extrapolations for odd and even 
lattice sites. The bars on the symbols give the "error" 
in the mean values, 5A(U;L). The statistical width in- 
troduced by the averaging procedure is small enough so 
that finite-size corrections are the dominant effect. 

The gap A([7), Fig. [j], is zero to within the error of 
the finite-size extrapolation (approximately the size of 
the points or less) for U < U c w W. As can be shown 
analytically in the large-t/ limit [Lql , the upper and lower 
Hubbard bands have width W for the case of a semiellip- 
tic density of states. In fact, above the transition we find 
A(L7 > 1.1IF) = U— W, with an accuracy of better than 
a few percent. An (almost) linear behavior for U > 1.2W 
was also reported in Refs. |K|[l4|]. Since 1 /[/-corrections 
do not contribute down to U = 1.1W, and finite-size cor- 
rections are largest below that value, we conjecture that 
A(U) =U-W might hold for U > U c = W. Allowing 
for the possibility of two critical interaction strengths, we 
give U c .i = (1.0 ± 0.1) W as a conservative estimate. 



In order to determine the quasiparticle weight Z(U), 
we need to resolve the momentum distribution in the 
vicinity of the Fermi energy. For even system sizes the 
closest point Ek to the Fermi energy E-p = is 5(^/2— ■ 
For odd system sizes there is one level right at the Fermi 
energy but the next one lies at £( 7r / 2 _2ir/L)- Therefore, 
we expect better results for Z(U) from an extrapolation 
of the even-size systems. In Fig. || we display the mo- 
mentum distribution n(ek) for two U values. The lines 
are fits to the Fermi-liquid form n(e < Q;U;L) = [1 + 
Z(U; L)]/2 + a(U; L){e/W) In (e/W) + f3{U; L)(e/W) + 
. . .] with n(—e) = 1— n(e) due to particle-hole symmetry. 
In Fig. H| we display Z(U;L), fitted by a second order 
polynomial in 1/L. The QJL) — > extrapolated value, 
Z(U), is displayed in Fig. |[ 
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FIG. 3. The energy-shell-averaged momentum distri- 
bution n(sk) for U = 5t, and U = St. The results for 
L = 6, 8, 10, 12, 14 are superimposed. 

We expect finite-size effects to be under good con- 
trol if, in the interacting density of states, we can de- 
tect the quasiparticle resonance, which, at half band- 
filling, represents ZL/2 electrons. To obtain a rather 
conservative bound we assume that we miss the reso- 
nance if it contains less than one electron in our largest 
system, Z mm = 2/L max = 1/7. If we further assume 
a reasonably regular behavior of Z(U) for U > 0.9 IF 
we find U c ,2 = (1-0 ± 0.1)W for the vanishing of the 
quasiparticle weight. This idea is supported by the 
fact that fourth-order perturbation theory, Z^(U) = 
1 - 1.308(C//IF) 2 + 0.6369(C//VF) 4 - . . . , quantitatively 
agrees with our results up to (7 w 0.6IF, above which 
higher-order corrections become important. 

This comparison is shown in Fig. |^, along with the 
results from iterated perturbation theory (IPT), the pro- 
jective self-consistent method (PSCM) and two 
exact diagonalization studies (EDI Q and ED2 Q) of 
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the self-consistency equation. Note that, aside from the 
standard fourth-order perturbation theory [PT 0(U 4 )], 
all of these results significantly deviate from the RDA 
results at U > W/2. 
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FIG. 4. The quasiparticle weight Z(U ; L) plotted versus 
1/L for various U values, fitted with a second-order polyno- 
mial in 1/L. 
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FIG. 5. The quasiparticle weight Z(U) within the RDA, 
fitted with a spline function and compared to various other 
approaches described in the text. Also included is U c ,2 ob- 
tained from the projective self-consistent method (PSCM). 

The IPT is a perturbative approach which may fail 
in the vicinity of the quantum phase transition n\. It 
is seen to be quantitatively unreliable even at moderate 
interaction strengths. The PSCM also predicts a con- 
tinuously vanishing quasiparticle weight and a jump dis- 
continuity in the gap. It requires the separation of high 
and low energy scales close to the transition. However, 
this assumption is in conflict with the mapping of the 



infinite-dimensional Hubbard model onto the dynamical 
mean- field equations |l(J . In the ED studies of the self- 
consistency equations, the energy levels of the effective 
medium are discretized, which limits the energy resolu- 
tion to Alu > 0.2W for n s < 10 bath levels Q. This 
energy resolution may be insufficient to distinguish the 
spectral weight of the quasiparticle resonance from spec- 
tral weight in the tails of the Hubbard bands accurately. 
The size of the quasiparticle weight and therefore the 
stability of the metal could thereby be overestimated. 

In this work we have applied the random dispersion 
approximation to calculate directly the gap and the mo- 
mentum distribution on finite systems. We have used 
conventional polynomial fits in 1 / L to extrapolate to the 
thermodynamic limit, and have found results in full ac- 
cordance with Mott's view of a continuous Mott-Hubbard 
transition with a single critical interaction strength. 
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